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Abstract. We have derived beryllium abundances in a wide sample of stars hosting planets, with spectral types 
in the range F7V— KOV, aimed at studying in detail the effects of the presence of planets on the structure and 
evolution of the associated stars. Predictions from current models are compared with the derived abundances 
and suggestions are provided to explain the observed inconsistencies. We show that while still not clear, the 
results suggest that theoretical models may have to be revised for stars with T e g <5500K. On the other hand, 
a comparison between planet host and non-planet host stars shows no clear difference between both populations. 
Although preliminary, this result favors a "primordial" origin for the metallicity "excess" observed for the planetary 
host stars. Under this assumption, i.e. that there would be no differences between stars with and without giant 
planets, the light element depletion pattern of our sample of stars may also be used to further investigate and 
constraint Li and Be depletion mechanisms. 
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£h ! 1- Introduction 



or Mayor & Santos 2001). One remarkable characteristic 



In the last few years we have witnessed a fantastic devel- 
opment in an "old" but until recently not very successful 
field of astrophysics: the search for extra-solar planets. 
Following the first success in exoplanet searches with the 
discovery of the planet around 51 Peg (Mayor & Queloz 
1995| ) , the number of known giant planets orbiting solar- 
type stars did not stop growing, being currently of 72 (in- 
cluding 7 multi-planetary systems)]^. 

Unexpectedly, the planets found to date do not have 
much in common with the ones in ou r own Solar S ystem 
(for a review see Marcy et al. 200C, Udry et al. 2001 



appears to be related with the parent stars themselves: 
stars with planetary companions are considerably metal- 
rich when compare d with single field dwa rfs (G onzalez 
19981; Santos et al. |2000t Gonzalez et al. bOOll; Santos 



et al. 2001a , |2001b ). To explain the observed difference 



Send offprint requests to: Nuno C. Santos, e-mail: 
Nuno . Santos@obs . unige . ch 

* Based on observations collected with the VLT/UT2 
Kueyen telescope (Paranal Observatory, ESO, Chile) using the 
UVES spectrograph (Observing runs 66.C-0116A and 66. D- 
0284 A), and with the William Herschel and Nordic Optical 
Telescopes, operated at the island of La Palma by the Isaac 
Newton Group and jointly by Denmark, Finland, Iceland, and 
Norway, respectively, in the Spanish Observatorio del Roque 
de los Muchachos of the Instituto de Astrofi'sica de Canarias. 

1 Counts of January 2002. For a list see e.g. table at 
obswww. unige . ch/~naef /who_discovered_that_planet .html. 



two main explanations have been suggested. The first and 
more "traditional" is based upon the fact that the more 
metals you have in the proto-planetary disk, the higher 
should be th e pro bability of forming a planet (see e.g. 
Pollack et al. 1996 for the traditional paradigm of plane- 
tary formation). Thus, in this case the "excess" of metal- 
licity is seen as primordial to the cloud that gave origin to 
the star/planet system. "Opposing" to this view, the high 
metal content observed for stars with planets has also been 
interpreted as a sign of the accretion of high-Z material 
by the star sometime after it reached the main-sequence 



(e.g. Gonzalez 1998; Laughlin 1996) 



Although recent results seem to favor the former sce- 
nario as the key process leading to the o bserved metal 
richness of stars wit h pla nets (Santos et al. 2001a , 2001b ; 
Pinsonneault et al. 2001 ), signs of accretion of planetary 
material have also been found for some planet hosts (e.g. 
Israelian et al. [2001a ; Laws & Gonzalez |2001 ). The ques- 
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Fig. 1. Spectra in the Ben line region (dots) for two of 
the objects observed, and three spectral synthesis with dif- 
ferent Be abundances, corresponding to the best fit (solid 
line) and to changes of ±0.2 dex, respectively. 



tion is then turned to know how frequent those phenomena 
happen, and to how much these could have affected the 
observed metal contents. 

One possible and interesting approach to this prob- 
lem may pass by the study of one particularly important 
element: beryllium (Be). Together with lithium (Li) and 
boron (B), Be is a very important tracer of the internal 
stellar structure and kinematics. Be is mainly produced 
by spallation reactions in the interstellar medium, whil e 
it is burned in the hot stellar furnaces (e.g. Reeves 1994 ). 
While most works of light element abundances are based 
on Li (the abundances of this element are easier to mea- 
sure), Be studies have one major advantage when com- 
pared with Li. Since it is burned at much higher temper- 
atures, Be is depleted at lower rates than Li, and thus we 
can expect to measure Be abundances in stars which have 
no detectable Li in their atmospheres (like intermediate- 
age late G- or K-type dwarfs). In fact, for about 50% of the 
known planet host stars no Li was detected (Israelian et 
al. 2001b). Furthermore, Li studies have shown the pres- 
ence of a significant scatter for late-type stars of similar 
temperature. This has also been observed in open clusters 
where all stars have the same age, and appears to be re- 
lated to the clusters's age, rotational velocities, pre-Main 
Sequence history, etc. (see, e.g., Garcia Lopez et al. 1994; 

1998| ; Jones et al. 1999), a fact that may 



Randich et al. 

complicate or even preclude a comparison. 

Given all these points, Be studies of planetary host 
stars can indeed be particularly important and telling. 
For example, if pollution has played some important role 
in determining the high-metal content of planet host stars, 



Fig. 2. Spectrum in the Bell line region (dots) for 14 Her 
(HD 145675), and two spectral synthesis with different Be 
abundances, corresponding to no Be (solid line) and to 
log7V(i?e)=0.5. Error-bars represent the photon noise er- 
ror. A conservative upper limit of 0.5 for the Be abundance 
was considered. 



we would expect to find a similar or even higher increase 
in the Be contents. This is basically due to the fact 
that planetary material is relatively poor in H a nd H e 
when compared to the star (e.g. Anders & Grevesse 1989 ). 
Furthermore, and unlike for iron, Be may be already a bit 
depleted in the stellar surface. Thus, the injection of plan- 
etary material into this latter could even be responsible for 
a more important abundance change in the Be abundance 
than in the iron content. If pollution has indeed played an 
important role, the net result of the fall of planetary mate- 
rial into the central "sun" would thus be that, for a given 
temperature interval, we should find that planet hosts are 
(in principle) more Be-rich that non-planet host stars. In 
other words, the analysis of Be abundances represents an 
independent way of testing the pollution scenario. 



Garcia Lopez & Peres de Taoro (199S) carried out the 
first Be measurements in stars hosting planets: 16CygA 
and B, and 55 Cnc, followed by Deliyannis et al. (2000). In 



order to continue to address this problem, we presentPj here 
a study of Be abundances in a set of 29 stars with planets, 
and a smaller set of 6 stars without known planetary com- 
panions. In Sects. H and ^ we present the observations and 
analysis of the data and in Sect.|| we discuss the results 
in the context of the planetary host stars chemical abun- 
dances, but also in terms of the Be depletion processes. 
We conclude in Sect.||. 



2 Pr elimin ary results were already presented in Garcia Lopez 
et al. (|200l|). 
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2. Observations and data reduction 

Observations of the Ben doublet at 3131 A were car- 
ried out during 5 different observing runs. Two of them 
made use of the UVES spectrograph, at the VLT/UT2 
Kucyen telescope (ESO, Chile). The obtained spectra 
have a resolution R= A/ A A of around 70 000 for one of 
the runs (66.C-0116A) and 50 000 for the other (66.D- 
0284 A), and the S/N ratios in the region around 3131 A 
varying between 30 and 250 (from now on we will re- 
fer to runs UVES (A) /(B), respectively). UVES spectra 
were complemented with data obtained using the UES 
spectrograph (R=55 000) at the 4.2-m William Herschel 
Telescope (WHT) and, in two different observing runs 
(A/B), with the IACUB echelle spectrograph (R=33 000) 
at the 2.6-m Nordic Optical Telescope (NOT), both at 
the Observatorio del Roque de los Muchachos (La Palma) . 
The observations, together with the corresponding observ- 
ing run and S /N ratios obtained are listed in Table |j} 

Data reduction was done using IRAF0 tools in the 
echelle package. Standard background correction, flat- 
field, and extraction procedures were used. For the 
UVES (A) run and for the UES data, the wavelength cali- 
bration was done using a ThAr lamp spectrum taken dur- 
ing the same night. For the UVES(B) and IACUB runs 
the wavelength calibration was done using photospheric 
lines in the region of interest. The final linear dispersions 
for the UVES(A), UVES(B), IACUB(A), IACUB(B), and 
UES spectra were 30, 17, 35, 35 and 17 mA pixel -1 , respec- 
tively, with rms values in the range 4xl0 -4 to 8xl0~ 3 A. 



3. Abundance Analysis 

3.1. Atmospheric parameters 

Garcia Lopez et al. ( 1995| ) carried out a detailed analysis 
of the sources of uncertainties regarding the determination 
of Be abundances. They concluded that the precision of 
the derived Be abundances is mostly dependent on the 
choice of the stellar atmospheric parameters. In particular, 
they are very sensitive to the adopted value for the surface 
gravity (log g). 

In order to limit the possible systematic errors in our 
determinations it is important, whenever possible, to use 
an uniform set of atmospheric parameters for all the pro- 
gramme stars. Wc thus decid ed to use the values listed by 
Santos et al. ( 2001a ; 2001b ), computed from an uniform 
and accurate spectroscopic analysis available for most of 
the stars studied in this paperQ. For three of the stars 
(HD870, LTD 1461, and HD 3823 for which no planetary 



3 IRAF is distributed by National Optical Astronomy 
Observatories, operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation, U.S.A. 

4 The list of atmospheric parameters in Santos et al. (2001a) 



companions were found to date), no parameters were avail- 
able, and we computed them using CORALIE or FEROS 
spectra, in the very same way as in Santos et al. (2001a). 
The values are listed in Table ^. As discussed by the au- 
thors, the errors in T e g are usually lower than 50 K, and 
errors in the microturbulence parameter are of the order of 
0.1 kms -1 . As for logg, the uncertainties are in the range 
of 0.10 to 0.15 dex[} 

3.2. Spectral synthesis 

The abundance analysis was done in standard Local 
Thermodynamic Equilibrium (LTE ) usin g a revised ver- 
sion of the cod e MO OG (Sneden |l975 ), and a grid of 
Kurucz et al. ( 1993 ) ATLAS9 atmospheres. Be abun- 
dances were derived by fitting synthetic spectra to the 
data, using the same line-list as in Garcia Lopez & Perez 
de Taoro ( |1998| ). While both Ben lines at 3130.420 and 
3131.065 A are present in our data, we only used the lat- 
ter, given the severe line blending in the region around 
3130.420 A (used only for checking the consistency of the 
fit). 

In the analysis, the overall mctallicity was scaled to 
the iron abundance. We then iterated by changing the Be 
abundance, the continuum placement and the Gaussian 
smoothing profile until the best fit for the whole spec- 
tral region was obtained (we fitted all the spectrum be- 
tween 3129.5 and 3132.0 A). When considered impor- 
tant (e.g. for tBoo), the smoothing function used was a 
combination of a Gaussian and a rotational profile; for 
these cases we used the v sini value determined from 
the width of th e CORAVEL cross-correlation dip (Bcnz 
& Mayor 1984 ). Three examples are shown in Figs. [I] 
and ||. The resulting abundances for all the objects ob- 
served are listed in Table ||. Here we use the notation 
logiV(.Be)=[Be]=log(Be/H)+12. 

3.3. Errors 

It is not simple to derive accurate uncertainties for mea- 



was determined by these authors for part o f the objects, or was 
t aken from the studies of G onzalez ( 199£ ) , Gonzalez & Laws 
(|l998|), and Gonzalez et al. (|l999|,|200l|) 



surements of Be abundances (Garcia Lopez et al. 1995). 
In this paper the errors were estimated as follows. We 
considered that from the errors of ±50 K in temperature 
and ±0.15dex in logg we can expect typical uncertainties 
around 0.03 and 0.06 dex, respectively. There are several 
OH lines blended with the Ben 3130.420 A line; changes 
in the oxygen abundance would also change the location of 
the pseudo-continuum in that region, affecting the overall 
fit. To take this into account, an error of 0.05 dex, associ- 
ated with the uncertainties in the oxygen abundances ex- 
pected for these stars, has been added. Other atmospheric 
parameters, like the metallicity [Fe/H] and the microtur- 
bulence, do not influence significantly the results, and we 
will conservatively consider that together they introduce 
an error of 0.05. Adding quadratically, these figures pro- 



5 Note that these small uncertainties are relative, and not 
absolute. 
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Table 1. Observing log for stars with and without planets. 
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V 
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May 
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HD 117176 
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May 
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HD 130322 
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8.0 
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14 Her 
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55 
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HD 195019 
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HD 10697 
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IACUB(B) 
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Oct. 


2001 


HD 12661 
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7.4 


IACUB(B) 


40 


Oct. 


2001 
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IACUB(B) 
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Oct. 


2001 
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v And 
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HD 120136 
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4.5 
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Aug. 
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HD 143761 
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5.4 
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120 


Aug. 
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51 Peg 
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100 


Aug. 


1998 


Stars without known planets: 












HD 870 
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7.2 


UVES(A) 


130 


Nov. 2000 


HD 1461 


HD 1461 


6.5 


UVES(A) 


120 


Nov. 2000 


HD1581 


HD1581 
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UVES(A) 


140 


Dec. 


2000 


HD3823 
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Oct. 
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HD 26965A 
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Oct. 


2001 
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7.2 


UVES(A) 


110 


Dec. 


2000 



duce an uncertainty of 0.09 dex, that was added to the 
error due to continuum placement and fit quality for each 
case, that we conservatively considered to be at least of 
0.10 dex. The final errors, listed in Table || together with 
the derived Be abundances, are of the order of 0.16 dex, 
and quite independent of the S/N of the spectrum. 

Note that we are interested in carrying out a differ- 
ential analysis, and thus the knowledge of the absolute 
temperatures and surface gravities is not very important. 
Rather, it is crucial that these values are all in the same 
"scale". 



4. Be in stars with planets 

There are many evidences indicating that the depletion of 
Be is connected with the rotational history of a star (e.g. 
Stephens et al. [1997 ). Although still not completely estab- 
lished, this history may be related to the presence or not 
of a (massive) proto-planetary disk (Edwards et al 
Stassun et al 



Strom 1994 



199S; Barnes et al. 2001: Rebull 



1993 



|2001| ; Hartmann pOOz] ). If true, this may result in differen ; 



depletion rates for stars which had different disk masses, 
and thus may or not have been able to form the now dis- 
covered giant planets. This could, in fact, have been the 
case for the pair of very similar dwarfs 16 Cyg A and B (the 
latter having a detected planetary companion), for which 
the Li abundances seem to be quite different (Cochran et 
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al. |!997t King et al. |1997| ; Gonzalez et al. |1998| ), but show 
similar amounts of Be (Garcia Lopez & Perez de Taoro 
199S). On the other hand, if pollution plays any role, we 



1.5 



might also expect to detect some differences between stars 
with planets and stars without planets concerning light el- 
ement abundances, and in particular concerning Be. 

Stars with planets might thus present statistically dif- 
ferent Be contents when compared with stars without (de- 
tected) planetary systems. 

In Fig. U we show a plot of Be abundances for the stars 
presented in this paper as a function of T c ff . Filled symbols 
represent stars with planets, while open symbols denote 
dwarfs with no planet companions found by radial- velocity 
surveys (i.e. measured in the context of the CORALIE sur- 
vey (Udry et al. 2000| ), and having no clear radial- velocity 
signature of a planet). We also included three objects from 



the study of Garcia Lopez & Perez de Taoro (1998) - trian- 



gles, namely 16CygA (HD 186408), 16CygB (HD 186427) 
and 55Cnc (HD 75732 A), for which the atmospheric pa- 
rameters have been computed using the same technique 
as for the stars presented here (Gonzalez 1998 ). In order 
to keep this work as an homogeneous comparative study, 
we have preferred not to introduce measurements for stars 
"wit hout" planets taken from other works (e.g. Stephens 
et al. 1997 ), derived using different sets of atmospheric pa- 
rameters and chemical analysis (e.g. equivalent widths). 

In the figure, stars with a surface gravity log g <4.1 
(probably slightly evolved) are represented by the 
squares. These include HD 10697, HD 38529, HD 117176, 
HD 143761, and HD 168443 (all planet hosts). Note also 
that HD 114762, a quite metal-poor dwarf ([Fe/H]=-0.6) 
with a brown-dwarf companion was also included amongst 
the planet hosts. 

Superposed with the measurements are a set of Yale 
beryllium depletion isochrones from Pinsonneault et al. 
( |l990|) for solar metallicity and an age of 1.7 Gyr. Given 
that it is not possible to know what was the initial Be 
abundance for each star, we have considered an initial 
log N (Be) = 1.28, i.e., b etween Solar (log N(Be)=1.15 - 
Chmielewski et al. |1975|) an d meteoritic (logiV(_Be)=1.42 
- Anders & GrevesseT98^) . 



Although still preliminary, a look at the figure shows 
that the current results argue against pollution as the key 
process leading to the met allicity excess of stars with plan- 
ets (s ee als o Santos et al. 2001a , 2001b and Pinsonneault 
et al. [2001D . Adding -50 earth masses of CI chondrites to 
the Sun, for example, would increase its iron abundance 
by about 0.25 dex (a value similar to the average differ- 
ence observed between stars with and without detected 
giant planets), and its Be abundance would increase by a 
slightly higher factor^. No such difference seems present in 
our data. In the same way, the results also do not support 



6 For this calculation we have used a Be/Fe ratio similar 
to CI chondrites, and cons idere d an iron content of ~20% by 
mass (Anders & Grevesse 198£ ) . Note also that the fall of a 
few jupiter-mass jovian planet (s) would not have a big effect 
on the iron or beryllium abundance of the star. 
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Fig. 3. Be abundances for stars with planets (filled sym- 
bols) and stars without detected planetary companions 
(open symbols) as a function of T c g . Circles represent the 
stars studied in this work, while triangles denote the three 
objects ta ken fr om the study of Garcia Lopez & Perez 
de Taoro ( [1998 ) . Stars with surface gravities log g <4.1 
are denoted by squares. Superposed are the Be depletion 
isochrones (case A) of Pinsonneault et al. ( 1990 ) for solar 
metallicity and an age of 1.7 Gyr. The Sun is represented 
by the Solar Symbol. From top to bottom, the lines repre- 
sent a standard model (solid line), and 4 models with dif- 
ferent initial angular momentum (tables 3 to 6). Typical 
error bars are represented in the lower-right corner of the 
figure. 



extra-mixing due to an eventual different angular momen- 
tum history of the two "populations" of stars. We note, 
however, that at this stage we are limited by the small 
number of comparison stars analyzed in this work. 

In this context, one particular star (HD 82943) de- 
serves a few comments. This planet host was re- 
cently found to have a near meteoritic 6 Li/ 7 Li ratio 



(Israelian et al. 2001a), better interpreted as a sign that 
planet/planetary material was engulfed by this dwarf 
sometime during its lifetime. Israelian et al. have sug- 
gested that the Li isotopic ratio of HD 82943 could be 
explained by the fall of either the equivalent of a 2 jupiter 
mass giant planet or a 3 earth mass terrestrial planet (or 
planetary material) . Considering that HD 82943 had ini- 
tially a solar Be abundance, the addition of this quan- 
tity of material would increase its Be abundance by about 
0.1 dex. A small "excess" of Be would thus be observed, 
but this value would be within the error bars of our mea- 
surements. As we can see from Fig.|^, the Be abundance 
of HD 82943 is not particularly high when compared to 
other stars in the plot (although it seems to occupy a po- 
sition near the upper envelope of the Be abundances). No 
further conclusions can be taken. 
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4.1. The log N(Be) vs. T cS trend 



e.g. 

basically no Be de- 



According to the standard models (without rotation 
Table 1 of Pinsonneault et al. |1990| ) 
pletion should occur for dwarfs in the temperature region 
plotted in the figure. This is clearly not compatible with 
the observations. However, it is already known that stan- 
dard models are not able to explain most of the observed 

19971 ) 



On 



behaviors of Li and Be (e.g. Stephens et al 
the other hand, models with rotation (Pinsonneault ct 
al. |lQQ0t )n 

were shown to be quite satisfac tory for stars 
with 6500K < T off <5600K (Stephens et al. |l997|) . But as 
we can see from the figure, they do not predict signif- 
icant burning at temperatures around 5200 K. We note 
that according to these isochrones, Be depletion is done 
at moderate rates once the star reaches the main-sequence 
(about 0.06 dex from 0.7 to 1.7 Gyr for a 0.9 M dwarf, 
and at a considerably lower r ate for higher ages following 
the Fig. 11 of Stephens et al. 1997 ), and thus we do not 
expect a crucial difference with the presented curves. 

In contrast with the models^, one very interesting de- 
tail is seen. If, for temperatures between about 5600 and 
6200 K, the dispersion in the log N(Be) is remarkably 
small, and the abundances seem to be consistent with the 
model predictions, the general trend for lower tempera- 
tures follows a slow decrease with decreasing temperature. 
Under the assumption that no difference exists between 
stars with and without giant planets, cither the trend is 
due to some metallicity or age effect, or it is simply telling 
us that the models are not able to reproduce the observa- 
tions for temperatures bellow ~5600K. This problem was 
also noted by Stephens et al. ( 1997 ), but no such extreme 
cases had been found, maybe because their "solar metal- 
licity sample" did not go down to temperatures lower than 
^5500 K. Note that the iron abundances for the stars in 
the plot are in the range from —0.6 to +0.5, and the vari- 
ation of the initial Be abundance with stellar metallicity 
seems to be q uite s mall in this interval (Rebolo et al. 1988 ; 
Molaro et al. 1997 ; Garcia Lopez et al., in preparation). 

Particularly interesting are the 4 cases for which only 
upper values for the Be abundances were found: HD 38529, 
HD 145675 (14 Her), HD13445 (G186), and HD 75732 
(55Cnc), all in the temperature interval between 5150K 
and 5700 K. A detail of the spectral synthesis for one of 
these stars can be seen in Fig.||. While for HD 38529 the 
Be depletion may be explained by the fact that this star is 
already leaving the main-sequence (its low surface gravity 
an d its p osition in the HR-diagramme - e.g. Gonzalez et 
al " 



2001 



show it to be evolved) for the three remaining 
objects no simple explanations seem to exist. 

One possible explanation for Be abundances of the 
most discrepant objects could be their ages. For exam- 



Improved Yale models (e.g. Deliyannis & Pinsonneault 



199S ) do not show an important difference. The 1990 mod- 
els were thus used in the current paper, since they are the only 
tabulated. 

8 We will assume for now that the initial Be abundance was 
similar (within the uncertainties) for all the stars in the plot. 
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Fig. 4. Lithium and Beryllium abundances for the stars 
studied in this paper. We note that for some of the 
stars, no Li abundances were found in the literature (see 
Table |). Symbols as m Fig. | Typical error bars are rep- 
resented in each panel (0.16 for Be and 0.12 for Li). 

pie, 55 Cnc and Gl 86 seem to be quite old, with isochrone 
ages in agreement with the oldest stars in the galaxy (ages 
were determined from theoretical isochrones of Schaller et 



al. |1992| , and Schaerer et al. |l993a| |l993bj ). This is not the 
case for 14 Her, with an age around ^6 Gyr. Furthermore, 
for 14 Her, there is a star (HD 222335, without detected 
giant planet) with about the same temperature, similar 
age (~5-6 Gyr) but much higher Be abundance. 

Can a metallicity effect explain the observed differ- 
ence? An increase in the metallicity (and thus in the 
opacities) has the effect of changing the convective en- 
velope depth of a star. It has been shown by Swenson et 
al. (1994) that even small changes in the oxygen abun- 
dances may have a strong effect on the Li depletion rates. 
Although it is unwise to extrapolate directly, a similar 
effect should exist concerning Be. In what concerns our 
case, most of the objects in the plot having temperatures 
lower than 5400 K are metal-rich with respect to the Sun. 
For example, 14 Her is one of the most metal rich stars 
in the sample ([Fe/H]=+0.50). But some counter exam- 
ples exist, like Gl 86, only about 100 K cooler, with a value 
of [Fe/H] =— 0.20. Unfortunately, no theoretical isochrones 
for Be depletion are available for these high metallicities, 
and it is not clear whether a difference of +0.5 dex in 
metallicity can or cannot induce significantly different Be 
depletion rates for this temperature range. On the other 
hand, if some unknown opacity effect could lead to system- 
atic errors as a function of temperatur e - e .g. the miss- 
ing UV opacity (Balachandran & Bell 1998 ) - we could 
maybe expect to measure some trend. However, we have 



N.C. Santos et al.: Beryllium abundances in stars hosting giant planets 



7 



no special reason to believe that this effect can lead to 
the observed trend. Furthermore, the missing UV opac- 
ity problem is far from being solved (if it exists at all); 
indeed, recent results seem to show that we can obtain a 
good fit of the solar spectrum without taking into account 
this extra-opacity (Allende Prieto & Lambert 200C). 

Although not conclusive, the evidences discussed above 
suggest that a bona-fide explanation for the observed Be 
abundances may pass either by some inconsistency in the 
models, or by some metallicity effect. But given the much 
higher number of planet hosts in the plot when compared 
with the non-planet hosts, we still cannot discard that the 
presence of a planet might be the responsible for the ob- 
served trendQ. Else, since Be depletion in the Pinsonneault 
et al. ( 1990[ ) models is strongly related to the angular mo- 
mentum lost, we would have, for example, to consider that 
there is some negative correlation between initial angu- 
lar momentum and/or angular momentum loss and stellar 
mass (which does not seem very plausible) . Unfortunately 
only the addition of more objects to the plot, and in par- 
ticular the determination of Be abundances for a larger 
sample of dwarfs with T e a <5400K will permit to better 
settle down this question. Such observations are currently 
in progress. 

It is also worth mentioning that on the other side of 
the T cff plane, both HD 120136 (tBoo) and BD -10 3166, 
positioned in the Li (and Be) dip region (e.g. Boesgaard & 



1.5 



King 2001), have a Be abundance that may be compatible 



with the models. 



4.2. Lithium vs. Beryllium 

Our data can also be used to further investigate the is- 
sue of Li and Be depletion in main sequence stars. Fig. [| 
compares Li and Be abundances for the same stars, and 
shows, as expected, that Li is burned much faster than 
Be, and its decline with temperature is much more clear. 
This is clearly expected from the models for a middle age 
solar type dwarf in this temperature regime. Only two 
points seem to come out of the main trend (in the Li 
panel). These are the cases of HD 10697 and HD 117176 
(70Vir). Th ese tw o stars deserve some particular discus- 
sion in Sect. [4. 2.1 . 



In Fig. H we plot Be vs. Li abundances for all the stars 
with T c ff <6300K for which both Li and Be abundances 
are available. The plot reveals no clear trends. Even 70 Vir 
and HD 10697 follow the main trend, probably attesting 
the normality of these two stars. Again, no significant dif- 
ference is seen between the two groups of stars (planet 
hosts and non-planet hosts). 

One interesting but expected detail in Fig. [5] is that 
there are no stars in the lower right part of the plot, a re- 
gion that would correspond to stars having depleted much 



Can the presence of a massive proto-planetary disk (able 
to form planets), or other variables connected to the process 
of planetary formation itself (e.g. timescale), induce different 
angular momentum histories for stars with different T e ff? 
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Fig. 5. Li vs. Be abundances for the stars in our sample 
having temperatures lower than 6300 K and for which both 
Li and Be abundances were available. Symbols as in Fig.^[ 
Typical error bars are represented in the lower-right corner 
of the figure. 



of their Be but that are still Li rich. This "gap" in the 
figure, plus the few stars in the upper right corner, sug- 
gests that there might be a correlation between Li and Be 



depletions (like the one found by Deliyannis et al. 1998 



Boesgaard et al. 2001 and more recently by Boesgaard & 
King [2001 for their sample of hotter dwarfs, but absent in 
the study of Dommgucz Herrera 1998, and in the sample 
of old cluster members of Randich et al. [2001 ). However, 
there is a significant number of stars in the upper left re- 
gion of the plot having Be determinations but only upper 
limits for the Li abundance (objects that are not present in 
Fig. 14 of Boesgaard et al.). This particular point in quite 
interesting, since it is, to our knowledge, the first time 
such a "population" of stars is discussed in the literature. 

In Fig.|| we have done the same plot as in Fig.||, but 
this time using different symbols as a function of the tem- 
perature of the object. As we can easily see from the plot, 
the objects in the three different temperature regimes cho- 
sen occupy clearly different positions in the diagramme. 
Stars with T e ff <5400 are positioned in the lower- left cor- 
ner. These correspond basically to objects having burned 
some Be, and that only have upper values for the Li abun- 
dance. In the upper-right corner of the figure are posi- 
tioned stars in the range T e ff>5900K. These correspond 
to objects that have essentially not burned any Be, having 
depleted only moderately their Li content. An intermedi- 
ate population, with 5400K<T e ff<5900K, having already 
considerably depleted its original Li but not its Be content 
is present in the mid- upper part of the plot. 
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three temperature intervals: T e g<5400K (filled circles), 
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Except for HD 38529, a sub-giant star (the object with 
the lowest Be abundance in the plot), all the other objects 
are well grouped according to their temperature range 
This is an interesting result, since it is giving us impor- 
tant clues about the temperature at which the onset of Li 
and Be depletion occurs in our metal-rich stars: if for Li 
depletion occurs already for temperatures around 5900 K 
for Be this only seems to happen at T e ff~5400K. Note 
however, that given the dispertion these limits have to be 
seen as approximate. 

Note also that we did not include t Boo in Figs.|| and 
|6|, since its the only object in the sample that is clearly in 
the Li-dip region. 

4.2.1. HD 10697 and 70Vir 

As mentioned above, HD 10697 and 70Vir do not fit 
into the mean trend in the lower panel of Fig.||. 
Considering as initial abundances of Li the "cosmic" value 
(log7V(£z)=3.31; Martin et al. |l994j ) and for Be the av- 
erage between meteoritic and solar (log N(Be)=l.28, as 
above), then the observed values for 70Vir indicate that 
Be (log N(Be)=0.82) is depleted about 3 times, while 
Li (logiV(Li)=1.76) is depleted by by a factor of 22. 
Nevertheless, such a pattern of Li and Be depletion is not 
so unexpected ( see e .g. Deliyannis & Pinsonneault 1997 ; 
Stephens et al. 1997 ), and is compatible with models of 
depletion including turbulent induced rotat ional mixing 
(see e.g. discussion in Stephens et al. 1997 ). But 70Vir 



confirmed by its low surface gravity logg=3.90), and it 
seems strange that an evolved star with its temperature 
may still have so much Li. 

One possibility would be to invoke Li and Be dredge up 
from a "buffer" below the former main-sequence convec- 
tive envelope (Deliyannis et al. 1990). But as discussed in 
e.g. Randich et al. ( 1999] ), this scenario does not seem to be 
supported by the observations. However, the fact that this 
star is evolved means that its current temperature is prob- 
ably different to the one it had during the main-sequence 
phase. With a near solar metallicity ([Fe/H]=— 0.03) and 
a mass of ~1.1 Mq, the temperature of 70 Vir would have 
been more close to 6000 K. A dwarf with this tempera- 
ture is not supposed to burn much of its Li (as we can see 
from the plot). So, we speculate that this star may have 
just left the main-sequence, and started to dilute (and/or 
burn) the Li (as well as Be) content present in its con- 
vective envelope. If the size of the convective envelope has 
increased "quite fast" , maybe it still did not have time to 
deplete/dilute all the Li, but i s already depleting part of 
its Be (Charbonnel et al. |2000D . 

The case of HD 10697 is in fact not very different from 
the one of 70 Vir. Both have the same mass (1.10 Mq), and 
similar surface gravities. The main difference resides in the 
fact that HD 10647 seems to have, within the errors, pre- 
served most of its Be content, while Li is depleted about 
23 times. But given the uncertainties in our Be abundance 
determination, and the unknown initial Li and Be content 
of the star, we cannot be certain whether the Be is really 
intact in the atmosphere of this object. The same argu- 
ments discussed for 70 Vir may thus be valid in this case. 

It is worth mentioning that another possibility to ex- 
plain the high Li abundances of these two stars would 
be to invoke a planet en gulfmen t (like in the case of 
HD 82943 - Israelian et al. 2001a ), that could simply re- 
sult from the migration of a former close-in planet due 
to tid al interactions with the evolved star (Rasio et al. 
1996 ). Unfortunately, no measurements of the Li isotopic 
ratio are avail able for 70 Vir. But recent models by Siess 
& Livio ( 1999] ) show that the absorption of a planet by a 
giant star would have the effect of increasing considerably 
its Li abundance, but only slightly the Be content, a situ- 
ation that could be compatible with the observed Li and 
Be abundances for this star. 

There is another interesting point about these two 
stars. Although they have very similar Li abundances, 
their Be contents differ by 0.56 dex. Given that both stars 
have similar metallicity (—0.03 and +0.16, for 70 Vir and 
HD 10697, respectively) their initial Be content is not ex- 
pected to be significantly different (Rebolo et al. |1988; 



is not a young object (it seems to be slightly evolved, as 



Molaro et al. 1997| ; Garcia Lopez et al., in preparation). 
However, we do not know up to which extent Be and 
Li abundances are uniform in the Inter-Stellar Medium. 
A different initial Li/Be ratio, together with the uncer- 
tainties in the Li and Be determinations, could probably 
lead to the observed effect. Else, we would need to invoke 
that the two stars have burned Li at the same rate, while 
the opposite happened for Be. But given that these two 
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stars are sub-giants, we cannot exclude that we are simply 
observing some evolutionary effect, eventually combined 
with slightly different initial Li and Be abundances. Note 
that the former object is a bit cooler, being positioned in 
the temperature region where the Be depletion seems to 
occur (see previous Sect.). 



5. Conclusions 

We have obtained Be abundances for a sample of planet 
hosts stars, and a smaller sample of stars without known 
giant planetary companions, aimed at studying the exis- 
tence of any significant difference between the two groups. 
Our data also allow us to further investigate Li and Be de- 
pletion among metal-rich stars. The objects have metal- 
licities between about —0.6 and +0.5 dex, and cover the 
temperature interval between ^5150 K and 6450 K. The 
main conclusions go as follows: 

— No particular difference was found between the Be 
abundances of stars with and without giant plane- 
tary companions in the temperature interval for which 
both groups were observed. This may be interpreted 
as an argument against the "pollution scenario" as 
the source of the high-[Fe/H] observed for planet host 
stars. However, given the low number of comparison 
stars available in this study, no strong constraints can 
still be set. 

— For stars with T e fi <5500K, the general trend of 
the Be abundances is not well described by models 
of light element de pletio n for solar metallicity (e.g. 

199C| ). These inconsistencies may 



Pinsonneault et al. 
be connected either to some metallicity effect, or to 
some problem with the models, but at the moment 
we cannot completely discard some relation with the 
presence of a planet. 

— For temperatures lower than 6300 K, no clear corre- 
lated Li and Be pattern is seen, contrary to what was 
found by other authors for hotter dwarfs (Deliyannis 
et al. |!998t Boesgaard et al. |200l[ ). 

— Our data permit to establish approximate temperature 
limits for the onset of strong Li and Be depletion for 
the metal-rich stars studied in this work. For stars with 
T e g <5900K, severe Li burning is seen, while for Be 
this limit is around 5400 K. 

A large survey for Be in planet and non-planet host 
stars is currently under way, and we expect it to give new 
insights into this timely and interesting subject. 
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5750 


4.10 


-0.29 





95±0.19 


UES 


1.30 


HD 145675 


5300 


4.27 


0.50 




<0.5 


IACUB(A) 


<0.7 


HD 168443 


5555 


4.10 


0.10 





97±0.15 


IACUB(A) 


<0.71 


HD 187123 


5830 


4.40 


0.16 





98±0.14 


IACUB(A) 


1.20 


HD 195019 


5830 


4.34 


0.09 


1 


27±0.15 


IACUB(A) 


<1.05 


HD 210277 


5575 


4.44 


0.23 


1. 


05±0.17 


UVES(A) 


<0.73 


HD 217014 


5795 


4.41 


0.21 





98±0.14 


UES 


1.30 


HD 217107 


5660 


4.42 


0.39 


1. 


01±0.17 


UVES(A) 


<0.86 


Stars from Garcia Lopez 


& Perez de Taoro ( 


1998 


): 






HD 75732 A 


5150 


4.15 


0.29 




<0.55 


- 


<0.04 


HD 186408 


5750 


4.20 


0.11 


1.10±0.17 




1.24 


HD 186427 


5700 


4.35 


0.06 


1.30±0.17 




<0.46 


Stars without known planets: 














HD 870 


5425 


4.59 


-0.03 


0.84±0.16 


UVES(A) 


<0.35 


HD 1461 


5785 


4.47 


0.18 


1.20±0.17 


UVES(A) 


<0.71 


HD1581 


5940 


4.44 


-0.15 


1.19±0.14 


UVES(A) 


2.35 


HD 3823 


5950 


4.12 


-0.27 


0.98±0.14 


UVES(A) 


2.44 


HD 26965A 


5185 


4.73 


-0.26 


0.56±0.20 


IACUB(B) 


<0.22 


HD 222335 


5310 


4.64 


-0.10 


0.87±0.18 


UVES(A) 


<0.35 



\ The values for the Li abundances were taken from Israelian et al. (2001b); we refer to this work for references. These authors 
used the very same set of stellar parameters as in this paper, and thus the listed Li and Be abundances are based on the same 
parameter scale. 

ft We refer to obswww.unige.ch/~naef/who_discovered_that_planet. html for the planet discovery references, 
ftt For this star (rBoo) the uncertainty in log N(Be) is higher due to the "large" v sini ~15ms~ . 



